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The Turbulent Destruction of Clouds - II. Mach Number 
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ABSTRACT 

The turbulent destruction of a cloud subject to the passage of an adiabatic shock is 
studied. We find large discrepan cies between the lifetime of the cloud and the analytical 
re sult offHartauist et al.l (|1986f ) . These differences appear to be due to the assumption 
in iHartQuist et al.l (|1986[ ) that mass-loss occurs largely as a result of lower pressure 
regions on the surface of the cloud away from the stagnation point, whereas in reality 
Kelvin-Helmholtz (KH) instabilities play a dominant role in the cloud destruction. 
We find that the true lifetime of the cloud (defined as when all of the material from 
the core of the cloud is well mixed with the intercloud material in the hydrodynamic 
cells) is about 6 x tKHD, where tKHD is the growth timescale for the most disruptive, 
long-wavelength, KH instabilities. These findings have wide implications for diffuse 
sources where there is transfer of material between hot and cool phases. 

The properties of the interaction as a function of Mach number and cloud density 
contrast are also studied. The interaction is milder at lower Mach numbers with the 
most marked differences occuring at low shock Mach numbers when the postshock gas 
is subsonic with respect to the cloud (i.e. M < 2.76). Material stripped off the cloud 
only forms a long "tail-like" feature if x ^ 10'^. 
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1 INTRODUCTION 

The interchange of material between dense cool phases 
and a hotter, more tenuous, external medium is a key 
phenomenon in astrophysics which influences the morphol- 
ogy and evolution of dilTuse sources over a wide range 
of spatial and energy scales. Clouds (also referred to as 
clumps, globules, and knots) may either accumulate ma- 
terial from, or lose material to, the surrounding medium. 
In the latter case, mass-loss can occur through hydrody- 
namic abl ation, or the rmal or photoionized evaporation 
(see, e.g., IPittardI |2007| . and references therein). The re- 
sponse of such multi-phase environments to the impact of 
winds and shocks is central to the investigation of feed- 
back in star and galaxy formation, and enc ompasses ob- 
jects such as tails behind mass-lo sing stars (e.g. Martin ot al.l 
I2OO7I ). planetary nebul ae fe.g. [ Matsuura et al . 2009), su- 
pernova remnants (e.g. Danforth. Blair fc Ravmondl I200ll : 



Strickland ot al.' '2000'; 'Cecil. Bland-Haw thorn" 
2002; Vcillcux, Cec il fc Bland-Hawthorn 20051') . and the in - 
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Levenson. Graham fc Waltersll2002l l. starburst regions (e.g. 



WestmoQuette et al.l 200?! . 20091 ) . starburst superwinds (e.g. 
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tracluster medium (e.g. lConselice. Gallagher fc Wvsell200ll ). 

Such interactions have been extensively investigated us- 
ing numerical hydrodynamics over the last two decades. A 
large literature of shock-cloud, wind-cloud, and jet-cloud in- 
teractions now exists. In it, the elTects of radiative cooling, 
thermal conduction, and ordered magnetic fields have all 
been considered. These studies have provided great insight 
into the behaviour and evolution of clouds subject to a va- 
riety of internal and external conditions. However, all suffer 
from a common and fundamental problem. In most astro- 
physical environments the interaction of shocks, winds, and 
jets with clouds takes place at high Reynolds numbers and 
is turbulent (see iPittard et al.ll2009l . for a review of the key 
physics). Unfortunately, the formation of fully developed 
turbulence is prevented by the artificial viscosity inherent 
in the hydrodynamical simulations. The turbulent mixing 
of cloud material into the surrounding medium is therefore 
limited, and the destruction of the cloud is hindered. 

The only tractable solution to this problem is to use a 
statistical approach to describe the turbulent flow. One pos- 
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sible method involves the use of a subgrid turbulent viscosity 
model, designed to calculate the properties of the turbulence 
and the resulting increase in the transport coefficients that 
the turbulence brings. In this way, the subgrid turbulence 
model emulate s a high Reyno l ds nu mber flow. An initial in- 
vestigation by iPittard et al.l (|2009l ) of an adiabatic, Mach 
10, shock-cloud interaction found increasing and significant 
differences between the results from a k-e turbulence model 
and an inviscid model as the density contrast between the 
cloud and inter-cloud medium, x, increased past 10'^. In ad- 
dition, the k-e turbulence model demonstrated better con- 
vergence in resolution tests, a feature which is particularly 
useful for simulations involving multiple clouds. The effect 
of a turbulent, as opposed to a laminar, post-shock medium 
sweeping over the cloud was also investigated. The strong 
"buffeting" that the cloud receives in this case results in its 
significantly quicker destruction. 

It is also surprising that previous shock-cloud studies 
have focussed almost exclusively on high Mach number in- 
teractions (see Table [ij). Where lower Mach numbers have 
been considered (e.g.T lFVagile et al] |2004| : iNakamura et al.l 
I2OO6I ). a detailed comparison to higher Mach number inter- 
actions is lacki ng. In this work we extend the shock-cloud 
investigation of lPittard et al.l (|2009l ) by examining the Mach 
number dependence of the interaction. We compare our nu- 
m erical results aga i nst th e analytical mass-loss rate formula 
of iHartguist et al] (|l986l l. and for the first time investigate 
in detail how well it describes the rate of destruction of the 
cloud. We also determine the conditions which are required 
for a long identifiable tail to be created behind the cloud. 

The structure of this paper is as follows. In Section [5] 
the numerical method and analysis are introduced. Section[3] 
notes the main stages in the destruction of a cloud by a high 
Mach number shock, and the Mach-scaling of the interaction 
is discussed. Our results are presented in Section |4l where 
a variety of simulations for strong and weak shocks are pre- 
sented and compared. Section [S] summarizes and concludes 
this work. 



2 THE NUMERICAL SETUP 

The calculations were performed in 2D axisymmetry using 
an Eulerian adapative mesh refinement (AMR) hydrody- 
namic code, with a linear Godunov s olver and piece-wise 
linear cell interpolation (see lFall3ll99lh . The entire compu- 
tational domain is covered by the two coarsest grids, G° and 
. The solution at each position is calculated on all grids 
that exist there, and the difference between these solutions is 
used to control refinement. 8 grid levels were used in total, 
with a factor of 2 refinement between each. In all simula- 
tions the coarse-grid resolution is equal to the radius of the 
cloud, with 128 cells per cloud radius on the finest grid. 
An advected scalar is used to distinguish between cloud and 
ambient material. 

The interaction is simplified by adopting a number of 
assumptions. The magnetic field is assumed to be too weak 
to be dynamically important. We also assume that radia- 
tive cooling can be ignored. This is valid if the cloud is 
small enough, and preserves the scale-free nature of the 
simulations. The efficiency of thermal conduction in mag- 
netized turbulent plasmas remains highly uncertain (see 



IPittard et al.|[2009l . and references therein), and for simplic- 
ity its effects are also ignored here. Finally, we also ignore 
self-gravity. It is unimportant in a cloud struck by a strong 
adiabatic shock (|Klein et al.lll993 ). but should be included 
in future work on mildly supersonic shock cloud interactions. 

The turbulence in the sub-grid k-e model is assumed 
to be fully developed, and has been calibrated by compar- 
ing the computed g rowth of shear layers with experiments 
(|Dash fc WoH|[l983l l. Turbulent energy {k, per unit mass) 
is generated by the action of the turbulent viscosity on the 
mean fiow and is converted to heat at the dissipation rate 
per unit mass, e. Since the turbulent energy resides mainly 
in large eddies, while the dissipation occurs in the small 
ones, one can think of k and e as describing the large-scale 
and small-scale turbulence respectively. Further details of 
the model implementation inclu ding the full set of eq uations 
that are solved can be found in IPittard et all (|2009h . 

The ratio of the turbulent and thermal energy densities, 
etb/eth, is ~ 10~^ in the initial post-shock fiow. Within the 
pre-shock medium (including the cloud), we set etb/cth = 
0.04. These values are low enough to not affect the initial 
dy namics of the i nterac tion, and are identical to those used 
bv IPittard et al.1 (|2009l ') in their investigation of a Mach 10 
shock-cloud interaction (where much higher levels of post- 
shock turbulence were also investigated). 

The numerical domain is set large enough so that the 
cloud is well dispersed and mixed into the post-shock flow 
before the shock reaches its edge. In this way various global 
quantities detailed below can be accurately computed. The 
shock propagates parallel to the axis of symmetry, which is 
the 2-axis. For clouds with a density contrast x ~ 10^ with 
respect to the ambient medium, the grid had an extent of 
< r < 24, with -910 < z 290 when M = 1.5, -510 < 
z 5? 20 when M = 3, and -510 < z < 6 when M = 10, where 
the unit of length corresponds to a cloud radius, rc. Smaller 
grids were used when x was lower. All calculations were 
performed for an ideal gas with 7 = 5/3, and were scaled 
so that the fluid variables have values reasonably close to 
unity. 

Clouds in the ISM do not have infinitel y sharp edges, 
so we adopt the density profile specified in IPittard et al.l 
(|2009l ): 



p{r) = Pamb[l/' + {l-^)r]], 



where 



V 



1 / Q - 1 
2(1 + ^ 



(1) 



(2) 



a = exp {min[20.0,pi((r/rc)^ — 1)]}, and r is the distance 
from the centre of the cloud. xf> is adjusted to obtain a specific 
density contrast for the centre of the cloud with respect to 
the ambient medium (x = Pmax/pamb). The parameter pi 
controls the steepness of the profile at the edge of the cloud. 
We set pi = 10 (i.e. a sharp-edged cloud), and place the 
cloud in pressure equilibrium with its surroundings. 

The evolution of the interaction is studied t h rough 
various integrated q uantities (se e Klein et al.l 1 19941 : 
INakamura et al.1 12OO6I : IPittard et al.1 l2009l ). These include 
the effective radii of the cloud normal to and along the axis 
of symmetry (a and c, respectively), its mass (m), its mean 
velocity ((uz), measured in the frame of the unshocked 
cloud), the velocity dispersions in the radial and axial 
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Table 1. A representative (but incomplete) summary of previous numerical investigations of shock-cloud and wind-cloud interactions. 
X is the density contrast of the cloud with respect to the ambient medium and M is the shock Mach number. The inte raction types are: 
SCS = s hock-cloud (single ); SCM = shock-cloud (m ul tiple); WCS = wind-cloud (single). T h e references are as fo ll ows: "ISto ne fc N ormanI 
lll992l'): iKiein et al.l jlQM): 'iMac Low et al.l ll 19941); 'iMellema. Kurk fc RottgerhiS ||2002|'); ^Fragile et Z l (1 20041'): ^ lOrlando et al.l ||2005^: 
^Nakamura et al.l ||2006|); van Loo et al.l ||2007|); lOrlando et al.l ||2008|'): - IShin et al.l ||2008|'I: *|Po ludnciiko ct al.l |20ol)~lGrego^^] 
l|200d l: " IPoludnenko et al.l l|2004l 'l: iMarcolini et al.l l|2005l ): feaga et al.l l|2007l l; ^Vieser fc Henslen t2007. ). 



Authors 


Interaction 
type 


Geometry 


Typical (max) 
resolution 


X 


M 


Cooling? 


Conduction? 


Magnetic 
fields? 


SN92'' 


SCS 


3D XYZ 


60 (60) 


10 


10 








KMC94'' 


SCS 


2D RZ 


120 (240) 


3,10,30, 
100,400 


10,100,1000 








MC94'^ 


SCS 


2D RZ,XY 


50 (240) 


10 


10,100 






✓ 


MKR02'* 


SCS 


2D RZ,XY 


200 (200) 


1000 


10 


✓ 






F04= 


SCS 


2D XY 


200 (200) 


1000 


5,10,20,40 


✓ 






005^ 


SCS 


2D RZ, 3D XYZ 


132 (132) 


10 


30,50 


✓ 


✓ 




N06S 


SCS 


2D RZ, 3D XYZ 


120 (960) 


10,100 


1.5,10,100,1000 








V07'' 


SCS 


2D RZ 


640 (640) 


45 


1.5,2.5,5 


✓ 




✓ 


O08* 


SCS 


2.5D XYZ 


132 (528) 


10 


50 


✓ 


✓ 


✓ 


SSSOS^ 


SCS 


3D XYZ 


68 (68) 


10 


10 






✓ 


PFB02'^ 


SCM 


2D XY 


32 (32) 


500 


10 








GOO' 


WCS 


3D XYZ 


16 (26) 


10,30,100 


1.5,3 






✓ 


PFM04'" 


WCS 


2D RZ 


128 (128) 


100 


10-200 


✓ 






M05" 


WCS 


2D RZ 


150 (150) 


100,500 


3,6.67 


✓ 


✓ 




R07° 


WCS 


3D XYZ 


77 (77) 


10 


242 


✓ 






VH07P 


WCS 


2D RZ 


33 (51) 


0.6,1.2, 
6.1 X lO'i 


0.3 


✓ 


✓ 





Notes: iPatnaude fc FesenI ||2005|') and [Cooler et al.l bOOgI) consider a cloud with substructure. Multiple-cloud interactions were also 

considered bv ll^agUe^t al. I l|2004l ) and iMeUolTet^ 1 l|2005l '). though this was not the main focus of the ir work. Simulations with 

contin u ous mass-injection (to m imic the des t ruction of very long l ived c louds) have been presented bv lFalle et al. I ||2002| ').I Pittard et al.l 
l|2005l '). lDvson et al.l l|2006l '). and lPope et al.l l|2008l ') . [Pittard et al.l l|2005l ') investigated the interaction of a wind with multiple 
mass-injection sources. 



Table 2. The dependences of the global cloud and core properties on the shock Mach number in subgrid turbulence calculations with 
cloud density contrasts of x = 10 (models "cl"), x = 10^ (models "c2"), and x = 10^ (models "c3"). The Mach number for each model 
is given by the number following the initial "m" in the model name - c alculations with Ma ch numbers of 1.5, 2, 3, 4, 6, 10, and 40 
were made. The time-dependent quantities are evaluated at t = tmix (cf. iPittard et ai]|2009l) . Values in parentheses are obtained from 
integrations over the "core" mass rather than the "cloud" mass. 



Model 


^dra 


y/tcc 


^mix/^cc 




a 


/rc 




c/ 


r-c 




c 


/a 




(P)/Pmax 




ml.5cl 


2 


97 


(3.04) 


(10.3) 


1 


90 


(1 


97) 


2.76 


[2 


74) 


1.45 


(1 


39) 


0.364 (0.543) 


0.292 


(0.305) 


m2cl 


1 


94 


(1.98) 


(7.95) 


1 


71 


(1 


78) 


2.68 


[2 


79) 


1.57 


(1 


57) 


0.452 (0.655) 


0.404 


(0.432) 


m3cl 


1 


40 


(1.47) 


(6.54) 


1 


52 


(1 


35) 


2.55 


[2 


75) 


1.68 


(2 


03) 


0.580 (0.909) 


0.471 


(0.472) 


m4cl 


1 


18 


(1.25) 


(6.11) 


1 


67 


(1 


58) 


2.46 


[2 


37) 


1.47 


(1 


50) 


0.649 (1.008) 


0.541 


(0.537) 


m6cl 


1 


04 


(1.13) 


(6.10) 


1 


77 


(1 


78) 


2.54 


[2 


70) 


1.44 


(1 


51) 


0.641 (0.883) 


0.604 


(0.629) 


mlOcl 


1 


04 


(1.11) 


(5.96) 


1 


66 


(1 


74) 


2.20 


[2 


40) 


1.33 


(1 


38) 


0.680 (1.030) 


0.579 


(0.629) 


m40cl 





91 


(1.01) 


(6.20) 


1 


78 


(1 


89) 


2.26 


[2 


42) 


1.27 


(1 


28) 


0.692 (1.044) 


0.599 


(0.657) 


ml.5c2 


6 


04 


(14.0) 


(8.41) 


2 


68 


(1 


46) 


9.63 


,1 


57) 


3.59 


(1 


07) 


0.046 (0.197) 


0.195 


(0.111) 


m2c2 


4 


43 


(4.70) 


(6.70) 


3 


11 


(2 


96) 


8.75 


;3 


35) 


2.81 


(1 


13) 


0.051 (0.109) 


0.354 


(0.307) 


m3c2 


3 


61 


(3.78) 


(6.04) 


4 


32 


(4 


64) 


6.90 


;i 


79) 


1.60 


(0 


38) 


0.056 (0.129) 


0.511 


(0.490) 


m4c2 


3 


44 


(3.49) 


(5.55) 


4 


09 


(4 


40) 


5.29 


[2 


20) 


1.29 


(0 


50) 


0.063 (0.121) 


0.536 


(0.535) 


m6c2 


3 


27 


(3.32) 


(4.98) 


3 


95 


(3 


28) 


4.62 


[1 


98) 


1.17 


(0 


60) 


0.073 (0.181) 


0.591 


(0.517) 


ml0c2 


3 


08 


(3.09) 


(4.86) 


3 


92 


(3 


35) 


4.12 


;3 


05) 


1.05 


(0 


91) 


0.075 (0.142) 


0.593 


(0.554) 


m40c2 


2 


96 


(3.01) 


(4.63) 


3 


88 


(3 


64) 


4.04 


[2 


29) 


1.04 


(0 


63) 


0.077 (0.163) 


0.586 


(0.542) 


ml.5c3 


9 


19 


(29.7) 


(13.3) 


5 


20 


(1 


66) 


65.5 


[A 


20) 


12.6 


(2 


53) 


0.0038 (0.029) 


0.199 


(0.073) 


m2c3 


8 


66 


(13.1) 


(10.5) 


3 


42 


(2 


55) 


57.4 


15 


86) 


16.8 


(2 


29) 


0.0051 (0.019) 


0.263 


(0.142) 


m3c3 


7 


92 


(9.73) 


(10.3) 


4 


31 


(4 


52) 


67.5 


19 


96) 


15.7 


(2 


21) 


0.0054 (0.012) 


0.430 


(0.360) 


m4c3 


6 


98 


(10.2) 


(8.34) 


3 


40 


(2 


30) 


61.2 


;6 


07) 


18.0 


(2 


64) 


0.0069 (0.027) 


0.394 


(0.213) 


m6c3 


6 


07 


(6.72) 


(6.95) 


4 


19 


(3 


78) 


42.1 


;5 


02) 


10.1 


(1 


33) 


0.0073 (0.019) 


0.449 


(0.333) 


ml0c3 


6 


84 


(7.15) 


(7.83) 


4 


01 


(2 


82) 


49.1 




88) 


12.2 


(2 


44) 


0.0078 (0.025) 


0.430 


(0.271) 


m40c3 


5 


29 


(5.63) 


(6.10) 


4 


47 


(4 


29) 


31.1 


15 


52) 


6.95 


(1 


29) 


0.0078 (0.018) 


0.467 


(0.366) 




1.34| 



21. 5| 
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Figure 1. Snapshots of the density distribution from a fc-e calculation of a Mach 1.5 adiabatic shock hitting a cloud with a density 
contrast of 10'^ with respect to the ambient medium (model ml.5c3). The evolution proceeds left to right and top to bottom with 
t = 0.02, 0.12, 0.27, 0.63, 1.34, 2.06, 4.22, 8.54, 14.3, and 21.5 tec- Note the different spatial scale used in the last three panels. The colour 
scale is logio p, and spans the range —4.5 (blue) to +1.0 (red). 

directions {5vi and Sv^, respectively), its volume {V), its 
mean density ((p)), and the total circulation produced 
(E). The whole of the cloud and the densest part of its 
core are distinguished by the v alue of the scalar v ariable k 
associated with the cloud (see iPittard et al.l |2009| ) . In this 
way, each global statistic can be computed for the region 
associated only with the core (identified with the subscript 



"core", e.g., ttcoro) or with the entire cloud (identified with 
the subscript "cloud", e.g., adoud)- 

The characteristic time for the cloud to be crushed 
by the shocks driven into it is the "cloud crushing" time, 
tec = K^^'^Tc/vt,, where v-t, is the velocity of the sh ock in 
the intercloud (ambient) medium (|Klein et al.lll993 ). Sev- 
eral other timescales are obtained from the simulations. The 
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Figure 3. Snapshots of the turbulent 

cvohition proceeds left to right with i 
routine. 



; energy per unit mass, k, from the k-e calculation with x = 10^ and M = 1.5 (model ml.5c3). The 
t = 0.12,0.27,0.63,1.34, and 2.06 tec- The white rectangular regions are artefacts of the plotting 
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time for the average velocity of the cloud relative to that of 
the postshock ambient flow to decrease by a factor of 1/e is 
defined as the "drag time", tdrag- The "mixing time", fmix, 
is defined as the time when the mass of the core of the cloud, 
m-corc, reaches half of its initial value. The zero-point of all 
time measurements occurs when the intercloud shock is level 
with the centre of the cloud. 

The diffusive nature of turbulence means that the cloud 
gradually disperses even when no shock is present. The time 
needed for the maximum density in the core of a cloud with 
X = 10 to drop to half its original value is ~ 150 tec, reducing 
to « 4 tec for a cloud with x = 10^- A comparison of these 
timescales to the cloud destruction timescales discussed in 
Section [4] reveals that the evolution of the cloud is always 
dominated by its interaction with the shock. 



3 STAGES AND "MACH SCALING" 

The main stages of an adiabatic, non-magnetized and non- 
conducting interaction of a high Mach number shock with 
a cloud h ave been described many times i n the literature 
(see, e.g., iKlein et all Il994l : IPittard et alll2009l 'l. Initially 
the cloud is compressed mainly in the 2-direction by the 
incident shock which propagates into the cloud and by a 
shock driven into the back of the cloud, and a bow shock 
propagates upstream into the ambient mediurrQ. The over- 



pressured cloud then expands downstream and laterally, and 
Rayleigh-Taylor (RT) and Kelvin-Helmholtz (KH) instabil- 
ities destroy the cloud and mix its material into the sur- 
rounding fiow. The cloud is more rigid, and is better able to 
resist the passage of the shock, as x increase^ 

The conditions behind a strong shock are virtually in- 
dependent of the sound speed ahead of the shock (e.g., for 
M = 10, the post-shock density is within a few percent of its 
value at M oo, while the normalized post-shock pressure 
is within a fraction of a percent). Since the inviscid, adi- 
abatic, non-magnetized and non-conducting hydrodynamic 
equations are invariant under the transformation 



V 



P 

M2' 



(3) 



with the position and density unchanged, the time evolution 
of the cloud is independent of the Mach number of the shock 
when expressed in units oit/tcc oc tM in the limit M — >■ oo. 
This is referred to as "Mach scaling" , and was d emonstrated 
to be reasonably valid for clouds with sharp (|Klein et al.l 
Il994h and with smooth l|Nakamura et al]|2006l ) boundaries. 

The sys tem of equations us ed for the fc-e calculations 
are shown in lPittard et al.l (|2009l ). They are invariant under 
the additional transformation 



M2' 



e 

M3' 



(4) 



so Mach scaling occurs in this case also when M ^ oo. 



The Mach numbers of the reflected and transmitted 
shocks, in the one-dime nsional planar limit, arc calculated in 
iMiesch fc Zweibe Unfortunately, in the strong, adiabatic 

shock limit, the analytical solution for the value of the Mach 
number of the reflected shock, M^, does not satisfy the original 



assumption that M'^ ^ 1, and so is known to be quantitatively 
incorrect. A 6*'^-order polynomial must instead be solved numer- 
ically. 

^ In fact, only clouds wi th x > 10"^ should be considered "rigid" 
jMiesch fc ZweibellfTool) . 
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At lower Mach numbers, the post-shock conditions are 
dependent on the shock Mach number and Mach scaling 
is not applicable. The evolution of the interaction is also 
different to the strong-shock shown below. 



4 RESULTS 

In this section we first examine the Mach number depen- 
dence of the interaction of a shock with a cloud of high 
density contrast (y = 10''). This extends the work in 
IPittard et all (|2009l l to a range of Mach numbers. The de- 
pendence of the interaction on x is then examined, after 
which we are able to draw some conclusions about the nec- 
essary conditions for material stripped from the cloud to 
form identifiable tails. We then examine three main themes: 
the Mach scaling of the interaction and the principal differ- 
ences which result when the post-shock flow is subsonic; the 
generation of turbulence in the interaction; and the mass- 
loss rate and the lifetime of the cloud. Some key parameters 
from the simulations are tabulated in Table [2] 



(measured in units of tec) as the Mach number increases 
(the M = 3 simulation forms a long well-defined tail by 
t — 8.54 tec)- This is because the stripping of material is 
more efficient at higher Mach numbers due to the faster 
growth of RT and KH instabilities. Fig. [2] also shows that 
the acceleration of the cloud increases with the Mach num- 
ber of the interaction, though at very high Mach numbers 
Mach scaling (Section [3]) holds. 

The development of turbulence in model ml.5c3 is 
shown in Fig.O where the turbulent energy per unit mass, k, 
is displayed, k is initially created in a thin turbulent bound- 
ary layer at the surface of the cloud where a region of high 
shear exists. The turbulent eddies are then advected by the 
fiow and a turbulent wake develops downstream of the cloud, 
with a setup time ~ tec- This beha viour is qualitatively sim- 
ilar to the Mach 10 case shown in lPittard erai] (|2009l ). al- 
though the peak turbulent intensity is smaller, and a su- 
personic vortex ring and its associated turbulence are not 
present. 



4.1 Cloud morphology and turbulence 

4- 1-1 Interactions with x = 10'^ 

Fig. [l] shows snapshots of the density distribution at dif- 
ferent times for the case of an adiabatic Mach 1.5 shock 
impacting a cloud with a density contrast x ~ 10^. The 
interaction is much milder than in the Mach 10 case (see 
IPittard et aLll2009^■ which resul ts in several major differ- 
ences (cf. Nakamura et al.ll2006l ): i) the postshock flow is 
subsonic with respect to the cloud, so a bowwave rather 
than a bowshock forms ahead of the cloud; ii) the compres- 
sion of the cloud is more isotropic; iii) a strong vortex ring is 
not produced; iv) the smaller velocity difference at the slip 
surface around the cloud limits the KH and RT instabilities 
which develop strongly in the M = 10 case; v) it takes much 
longer for the cloud to be mixed into the surrounding flow 
and for it to accelerate to the intercloud postshock speed. A 
comparison with the results from an inviscid code (without 
the subgrid turbulence model) reveals that the evolution is 
overall very similar. T his was also th e case for the Mach 10 
interaction studied in IPittard et al] 12009"), though signif- 
icant differences occur if the post-shock medium sweeping 
over the cloud is very turbulent. 

Direct comparison of the density structure at f = 
5.66 tec for interactions with Mach 1.5, 2, 3, 6 and 10 shocks 
are made in Fig. [2] There are significant differences between 
the density distributions from the low Mach number inter- 
actions on the one hand, and the high Mach number inter- 
actions on the other. At high Mach numbers [M > 2.76) 
the post-shock flow is supersonic, and a bowshock devel- 
ops around the cloud. The shock driven into the cloud is 
stronger, and less symmetric. The faster flow speed past the 
cloud results in faster growth of RT and KH instabilities, 
and more rapid acceleration of the cloud. 

The M = Q and M = 10 simulations show dense elon- 
gated tails at t — 5.66 tec, whereas the lower Mach number 
simulations do not. Since long tails do eventually form even 
in the Mach 1.5 interaction, it is apparent that material 
stripped off the cloud forms a longer tail at a given time 



4-1-2 Interactions with x = lO'^ and x = 10 

Clouds with density contrast x ~ 10^ are less rigid obstacles 
to the passage of a shock than clouds with x ~ 10'^. This 
can easily be discerned by comparing snapshots of the mass 
density at t = 4.85 tec in interactions of shocks of various 
Mach numbers with clouds of density contrast x ~ 10^ as 
shown in Fig. [J] to the snapshots in Fig. [2] of interactions 
with X = 10^ 

In fact, iKlein et al. noted that the nature of the 

interaction appears to change with x'- a^t low x the cloud is 
liable to break up into several large fragments, whereas at 
large x the core remains more intact and material is stripped 
from its surface (this is a less effective mixing mechanism). 
This observati on is partially s uppor ted by the density snap- 
shots shown in lPittard et al.l (|2009l ). where the last panel in 
each of Figs. 9, 10 and 4 shows snapshots at similar times 
for a Mach 10 shock interacting with a cloud where x ~ 10, 
lO'^, and 10'^, respectively. However, a closer examination 
of the statistics of the interaction reveals that in almost all 
of the cases investigated in this work, the x = 10^ cloud 
is destroyed faster than clouds with x = 10 and 10^ (the 
exception is the M = 1.5 simulation, where the x = 10 
cloud is actually destroyed faster). The generally rapid de- 
struction of clouds with X = lO'^ appears to be due to the 
fact that such clouds undergo a more rapid increase in their 
transverse radius as they expand following the initial com- 
pression caused by the passage of the shock (see Fig. [5]). 
This makes the cloud more vulnerable to fragmentation. 

The interaction of an adiabatic shock with a cloud of 
density contrast x = 10 has been discussed many times in 
the literature, which in the inte r est of brevity we do not 
repeat here. However, iKlein et al.l (|l994l ) notes that the axial 
stretching of the cloud increases with x- This is important 
insofar as we do not see a well-defined tail in any of our 
simulations when x ^ lO'^- It appears, therefore, that well- 
defined tails only form in interactions with density contrasts 
X ^ 10^ , though in such cases they can form across a wide 
range of Mach numbers. 
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Figure 5. Time evolution of the effective transverse radius of the cloud, acloudi for various Mach numbers and density contrasts: (a) 
X = 10, (b) x = 102, (c) x = 10-^ 
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Figure 6. Time evolution of the volume of the core of the cloud, Vcorc, normalized to its initial value, for various Mach numbers and 
density contrasts: (a) x = 10, (b) x = 10^, (c) x = 10^- 



4.2 Mach scaling and principle changes for 
subsonic post-shock flow 

In this section we further examine how the interaction 
changes with the Mach number, and the major differences 
which occur when the post-shock flow is subsonic. 

4-2.1 Cloud shape, volume, and mean density 

Fig- El shows that the initial compression of the cloud in the 
transverse direction is weakest in the M = 1.5 simulation, 
but comparable in the M = 3 and M = 10 simulations. 
The transverse expansion which follows as the cloud seeks 
to re-establish pressure equilibrium is typically slower in the 
M = 1.5 simulations, but occurs at a similar rate in the 
M = 3 and M = 10 simulations. 

That the initial compression of the cloud is weaker at 
low Mach numbers can also be discerned from Fig.[6lja). For 
low density contrasts (x 10) the core of the cloud main- 
tains a larger volume than when hit by a stronger shock. 
However, clouds with higher density contrast hit by a strong 
shock undergo a violent re-expansion after their initial com- 
pression, causing their core volume to exceed their counter- 
parts hit by weaker shocks (see panels (b) and (c) in Fig. [6]). 
At later times the volume of the cloud core decreases as 
the core material is gradually ablated and mixed into the 
surrounding flow. 

The reduced compressions behind low Mach num- 
ber shocks relative to their stronger counterparts cause 
lower peak densities of the clouds and their core mate- 
rial, (pcioud>//5max and (pcoro ) /pmax , respectively (see Fig.[7|. 
And although the decline in these quantities over the period 
i i$ 5 fee is slower in the Mach 1.5 simulations, due to the 



slower stripping of material into the flow, the flnal value of 
(Pcioud)/pmax that is reached is lower in the Mach 1.5 sim- 
ulations than in the higher Mach number simulations, since 
the density of the postshock ambient material is lower in 
this case. This is particularly noticeable in Fig. Ha). 

4-2.2 Cloud velocity 

At lower Mach numbers there is a gentler acceleration of the 
cloud and its core in the axial direction, shown by the mean 
velocity statistics, (dz, cloud) and (wz.coro), respectively (pan- 
els (a-c) and (d-e) in Fig. |H} . This reflects both the weaker 
shock that is initially driven into the cloud, and the slower 
speed of the post-shock flow past the cloud. The acceleration 
of the cloud to a velocity of 1/e times the ambient postshock 
flow speed is discussed further in Section [4.2.61 

4-2-3 Cloud mass 

The slower post-shock flow in the simulations with lower 
Mach numbers affects the growth rate of RT and KH in- 
stabilities. This in turn affects the speed at which mate- 
rial is ripped from the surface of the cloud and the time 
for the cloud to mix into the flow. Hence, mcorc declines 
more rapidly in simulations with higher Mach numbers (see 
Fig. [9]). In high Mach number interactions the core com- 
pletely disappears shortly after fmix, whereas at lower Mach 
numbers material may be identified as originating from the 
core for times considerably after tmix. Also of note is that 
the slope of rricore in the M = 1.5 simulations appears to get 
shallower with increasing x (Fig.[9ll. The maximum slopes of 
the M = 3 and M = 10 models are particularly steep when 
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Figure 7. Top panels: Time evolution of the mean density of the cloud, (pcloud)i normalized to the initial maximum cloud density, for 
various Mach numbers and density contrasts: (a) x = lOi (b) X = 10^, (c) x = 10'^- Bottom panels: As top, but for the mean density of 
the core, (pcorc)- 
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Figure 8. Top panels: Time evolution of the cloud mean velocity, (i'z)cloudi for various Mach numbers and density contrasts: (a) x = 10> 
(b) X = 10^1 (c) X = 10'^- Bottom panels: as the top panels, but for the core mean velocity, (?^z)core- 



X ~ 10^, which is hkely related to the rapid transverse ex- 
pansion of the cloud radius, and subsequent strong mixing, 
taking place in these models (see Section |47L2]) . 

The evolution of the cloud and the way that its mate- 
rial mixes with the surrounding gas can also be studied via 
its mass distribution in density or velocity space. Figs. 1101 
and II II show this evolution for clouds with x = 10^. Fig. 1101 
shows the fraction of cloud mass over a range of density 
bins of width 0.1 pmax. The first subfigure is shown for 
t = 1.34 tec, when the transmitted shock has swept through 
the entire cloud (see, e.g.. Fig. [l}. In the M — 1.5 simu- 



lation, the weak shock transmitted into the cloud initially 
compresses the gas by a factor of 1.7. Slightly higher val- 
ues are subsequently obtained as the shock converges on the 
cloud centre. In contrast, the stronger shocks in the M = 3 
and M = 10 simulations produce highe r compressions whic h 
can exceed a factor of 10 (see Fig. 19 in lPittard et al.ll2009l '). 
At later times the density of the cloud material drops slowly 
towards the ambient post-shock density as the cloud contin- 
ues to expand into and mix with the surrounding gas. 

Fig. Uni also reveals that flatter mass distributions are 
obtained during higher Mach number interactions, whereas 
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Figure 9. Time evolution of the core mass, mcoro, for various Mach numbers and density contrasts; (a) x = lOi (b) X = 10^; (c) X = 10'^ 
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Figure 10. Time evolution of the mass density distributions for models with x = 10^ and varying shock Mach numbers. The histograms 
indicate the cloud mass contained within each density bin of width O.lpmax, normalized to the total cloud mass, nric. From left to right 
the mass fractions were computed at t = 1.34, 2.78 and 5.66 tec- 
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Figure 11. Mass distributions as a function of (normalized to the post-shock velocity in the ambient medium, u^^^) for models with 
X = 10"' and varying shock Mach numbers at t = 5.66, 14.3 and 21.4 tec The histograms denote the mass contained within a velocity 
bin of width 0.01 «ics for a line of sight parallel to the z-axis. The integrated mass is the mass of the cloud, rric- 




Figure 12. Time evolution of the circulation, Fcirci for various Mach numbers and density contrasts: (a) x = lOi (t>) X = lO^i (o) 
X = 103. 
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the distribution is more concentrated and peaky in simu- 
lations with lower Mach numbers. Two distinct clusters of 
mass at densities above and below 2 pmax can be seen in the 
M = 1.5 simulation a.t t — 2.78 tec- This is not just a conse- 
quence of the higher Mach number cases being more evolved 
- the mass distribution in the M = 1.5 simulation never ap- 
proaches the flatness seen in the higher Mach number cases, 
even at much later times. 

The flatter mass distribution in the high Mach num- 
ber cases reflects differences in the dynamical state of the 
cloud core with the Mach number of the interaction. At 
t ~ 1.34 tec, the cloud hit by t he Mach 10 shock is vio- 
lently re-expanding (see Fig. 4 of iPittard et al.|[2009l '). This 
expansion is supersonic, which leads to a relatively low den- 
sity in the central region of the cloud, with higher densities 
in the shocked region around the periphery of the cloud. At 
t « 2.78 tec, the cloud core is relatively dispersed, with fin- 
gers of somewhat lower density material around its edges. 
Therefore, at both of these instances, the density of mate- 
rial in the core of the cloud spans a wide range. In contrast, 
in the Mach 1.5 case. Fig. [T] shows that at t = 1.34 tec, the 
core density of the cloud is approximately constant. Since 
the interaction is much milder, the gentler rebounding of 
the cloud in the Mach 1.5 case does not lead to such a wide 
range of densities, and therefore a more concentrated and 
peaky distribution is obtained as seen in Fig. 1101 

Fig. [TT] shows the mass distribution function in velocity 
space within the entire cloud integrated along the z-axis. 
The histograms indicate the fraction of the cloud mass con- 
tained within a velocity bin of width 0.01 Uics, where itics is 
the post-shock velocity of the ambient medium. Initially the 
cloud is at rest, and its velocity is zero. When the material 
within the cloud is fully mixed into the surrounding flow its 
velocity is equal to Uics (a small fraction of mass exceeds this 
value due to turbulent motions). The slower acceleration of 
material in the AI — 1.5 simulation is again clearly evident. 



Table 3. Analytical estimates for the total circulation and its 
various components produced by the postshock flow. In each case 
the calculation of the value uses the drag-time for the core as given 
in Table [H THng = 0.417,0.667, and 0.743 when M = 1.5,3 and 
10, respectively. In the fifth column the numerically determined 
peak values of the total circulation, — rtot)p, are given, while the 
sixth column notes the dominant componont(s). 



Model 


Rshock 


Tpost 


— Ttot 


— Ttot )p 


Dom. comp. 


ml.5cl 


0.85 


0.42 


0.84 


1.60 


Tshock and F, 


m3cl 


1.37 


0.52 


1.22 


1.95 


Tshock and F, 


mlOcl 


1.52 


0.48 


1.26 


1.75 


I^shock 


ml.5c2 


1.12 


6.06 


6.76 


6.76 


Fpost 


m3c2 


1.80 


4.20 


5.33 


7.48 


Fpost 


ml0c2 


2.00 


4.26 


5.52 


6.60* 


Fpost 


ml.5c3 


1.20 


40.7 


41.5 


37.2 


Fpost 


m3c3 


1.94 


34.2 


35.4 


36.5 


Tpost 


ml0c3 


2.15 


31.2 


32.6 


34.2 


Fpost 



Note: (*) The absolute value of the total circulation in model 
ml0c2 was still rising at the time the simulation was stopped. 

tal circulation can be seen in Fig. I12r a'). The passage of 
the shock over the cloud causes the initial rise to maxi- 
mum (Fshock), with the formation of the vortex ring be- 
hind the cloud (Fring) causing the subsequent drop. The 
post-shock flow past the cloud generates further vorticity 
which increases the circulation after this minimum (Fpost)- 
Figs. I12r b) and (c) demonstrate the increasing dominance 
of Fpost with increasing x- 

The predicted total circulation agrees better with the 
numerical results when the core (rather th an the cloud) 
value is used for tdrag in Eq.[6](as also found bv lPittard et al.l 
l2009t ). Table O shows the predicted and measured values. 



4-2.4 Circulation 

A central aspect of the interaction of a shock with a cloud 
is the developme nt of vorticity, lj, and circulation, F (see 
iKlein et al.lll99j ). Production of the latter can be divided 
into three main components: at the interface between the 
cloud and the surrounding flow by the initial passage of the 
shock (Fshock) and the subsequent postshock flow (Fpost), 
and at the triple points associated with the Mach-reflected 
shocks behind the cloud (F^ing). For the general case of a 
shock of Mach number M and a ratio of specific heats 7, 
these components are: 



6 / Af2 



7-f 1 y IVP 
1 ( hp 



1 - 



post 



(7 + 1)^ 
2 / Af 2 



A/2 
1' 



A/2 



Vbrc, 



(5) 



^CC 



7 + 1 \^ A/2 



"^^^ > ^br-c, (6) 
(7) 



Note that the circulation generated by the vortex ring is 
positive, while the passage of the shock and the postshock 
flow generates negative circulation. 

The contribution of these three components to the to- 



4-2.5 Energy evolution 

Figs. [13] and [14] show the gain in kinetic and thermal en- 
ergy of the cloud material, and the growth (and subsequent 
decay) of turbulent energy for simulations as a function of 
Mach number and cloud density contrast. The cloud mate- 
rial should eventually acquire the same kinetic and thermal 
energy density as the ambient medium, and the turbulent 
energy should dissipate as heat. The ratio of kinetic to ther- 
mal energy in the postshock flow of an adiabatic shock is 
/5k//5th = 7A/ps/3, where Afps is the postshock Mach num- 
ber measured in the frame of the undisturbed, upstream, 
ambient medium. It can be easily shown that 



A/ps = 



2 (A/2 



^[27A/2-(7-l)][(7-l)A/2 + 2] 



(8) 



and A/ = 10 shocks, at late times we expect that _Ek/£th = 
0.145, 0.607, and 0.954 respectively. Note that E^^/Eth = 1.0 
in the strong shock {M — ^ 00) limit. 

The plots of the kinetic energies of the cloud material 
in Fig. [13] are normalized against the values which the cloud 
material should eventually obtain once it fully mixes into, 
and becomes indistinguishable from, the post-shock flow. 
As can be seen, the timescales for this can be very long. 
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especially when the Mach number of the shock is low and the 
cloud density contrast is high (when M = 1.5 and x = 10'', 
the cloud material has attained less than 60 per cent of its 
fully-mixed kinetic energy at t = 25 tec)- 

Fig. [13] shows that the cloud's rate of gain of kinetic 
energy at a given value of x is similar in the M = 3 and 
M = 10 models, and in comparison is considerably slower 
in the M = 1.5 models. The same behaviour is seen for the 
thermal energy. In shock-cloud interactions with M — 1.5 
the rate of energy gain is very gradual and almost linear in 
time. In contrast, when M ^ 3, there is rapid transfer of 
kinetic energy to clouds with x ~ 10, with a much slower 
transfer of thermal energy such that Ey^/Eth can be signifi- 
cantly greater than unity over a substantial time period, be- 
fore finally dropping back below unity as it approaches its 
asymptotic limit. The transfer of kinetic energy to clouds 
with X = 10^ fi-nd 10'^ shows a noticeable and relatively 
short-lived increase in its rate when t ~ tdrag for the cloud 
(this is also when t ~ fmix), after which it slows. Fig. [14] 
shows that the cloud's rate of gain of kinetic and thermal 
energy at a given Mach number is initially much more rapid 
for lower values of x- 

The evolution of the turbulent energy is discussed in 
Section 14.31 



4.2.6 

Fig. [15] shows the Mach number and x dependence of tdrag 
and tmix, two of the important timescales which characterize 
the evolution and destruction of the cloud. The results of 
least-squares fits to these plots are noted in Tables [4] and [5] 
Both tdrag and tmix are relatively constant for a given x '<^^ 
Mach numbers above 4. This refiects the onset of the "Mach- 
scaling" mentioned in Section [3] In contrast, fdrag and tmix 
both rise sharply at lower Mach numbers: cloud drag and the 
mixing of material from the cloud are less efficient with low 
shock Mach numbers. This is because: i) the density jump 
across the shock and the speed of the postshock intercloud 
fiow past the cloud are lower when M is smaller; and ii) the 
growth rate of KH instabilities is slower. 

The cloud is accelerated by two processes. First, the 
shock driven into the cloud accelerates it to a speed Vs. Fur- 
ther acceleration then occurs as the shocked intercloud gas 
fiows past the cloud, until they have the same velocity. For 
large x l^i^ second stage dominates and Vs/vh is small. Solv- 
in g the equation of motion for the cloud (see, e.g., Eq. 2.5 
in lKlein et al.|[T994l ) in the strong shock limit for 7 = 5/3 
gives 



drag(s) 



1.53 



X 



-tc 



(9) 



where Cp is t he drag coefficient (see Section 2.1 of 
iKlein et al.lll994l ). In the general case we find 



*drag(g) — 2.3 



(7 - 1)M2 + 2 X 



1/2 



" ice • 



(10) 



(Af2 - 1) Cd 

For both Egs. l^and llOl the cross-sectional area of the cloud 
is assumed to remain constant. This is a poor approximation 
which leads to an overestimate of the drag time, since the 
shock causes the cloud to rapidly expand in the transverse 
direction (Fig. [5]). Nonetheless, it is useful to compare the 
ratio tdrag(g) /idrag(s) wlth M. For 7 = 5/3 and M = 1.5, we 



find that tdrag(g) /idrag(s) ~ 4. Comparing to Fig. [TSt a') we 
see that the ratio of tdrag when M = 1.5 to when M = 40 
is 3.3, 2.0, and 1.7 for x = 10, 10^ and 10^ respectively. 
Hence the Mach number dependence of the drag time scales 
nearly as expected when x = 10. However, the drag time 
for low values of M is increasingly overestimated compared 
to the simulations as x increases. This divergence is at least 
partly due to differences in the expansion behaviour of the 
cloud: comparing simulations ml.5c3 and ml0c3, one sees 
that Scioud is significantly greater in the former until t ~ 7tcc 
(see Fig. [5]:). 

At face value the fits in Fig. [I5ja) to the numerically 
obtained drag times imply that Co ~ 5, 5, and 10 for x = 10, 
10^, and lO"' respectively. These values are much larger than 
one expects (the true va lues of Co are likely to be < 1 - see 
iLandau fc Lifshit3ll959l . for solid bodies). This difference is 
again caused by the fact that the cross-sectional area of the 
cloud does not r emain constant. 

Prior to the iKlein et al.l (| 19941 ) paper, the naive expec- 
tation was that the mixing time of the cloud would be com- 
parable to the time taken for the cloud to sweep up a column 
of postshock mate rial of similar mass to the cloud. However, 
iKlein et al.l l| 19941 ) discovered that this was not consistent 
with their results, since there was no evidence for the mix- 
ing time scaling as x^^^^-cc, and proposed instead that the 
relevant time was the timescale for KH instabilities to frag- 
ment the cloud. The KH growth time is 



^KH 



(-) F^^^ 



(11) 



where kx is the wavenumber of the perturbation. Longer 
wavelengths {kxVc ~ 1) are the most disruptive. Assum- 
ing that the relative speed between the postshock fiow 
and the cloud is the post-shock flow speed (i.e. that the 
cloud is initially stationary after passage of the shock), then 
t'rci/'^b = (1 — po/pa), whcrc pQ aud ps are the preshock and 
postshock densities of the ambient medium, respectively. 
Setting kxVc — 1, one then obtains 



(7 + l)M% 
'™ ~ (2M2-2)*^ 



(12) 



While setting tmix = ^kh does not give the correct magni- 
tude for tmix (cf. Fig. I15p . the relative change with M (a 
factor of 1.8 between M — 40 and M — 1.5) is close to what 
is observed in our numerical model. 

Fig.[T5]also shows that the mixing time of the core, tmix, 
is always greater than the drag time of the cloud, tdrag- The 
ratio of fmix/tdrag is also dependent on this ratio being 
greatest at low values of x, but declining with increasing x 
as the relative efficiency of mixing relative to acceleration 
increases. For clouds with x ~ 10, the numerical simula- 
tions reveal that fmix/idrag increases from 1.39 at M = 1.5, 
peaking at 1.67 at M = 3, and thereafter declines to 1.56 at 
M = 40. When x = 10^, tmix/tdrag declines from a value of 
1.45 at M = 1.5, to 1.15 at M = 40. 

The values of fdrag and tmix display greater scatter 
about the least squares fit when x = 10'' compared to the 
fits at lower values of x- This is due to large scale RT insta- 
bilites in these models which randomly and spontaneously 
fragment the cloud. The values of tdrag and tmix are generally 
higher in inviscid calculations, and display larger scatter. 
The extra viscosity that the subgrid turbulence model im- 
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Figure 13. The evolution of the cloud kinetic energy (normalized 
to its asymptotic value once the cloud is fully mixed into the post- 
shock flow) in models with x = 10 (top), 10^ (middle), and 10^ 
(bottom). The ratio of the cloud's kinetic to thermal energy and 
its turbulent energy at subgrid scales as a fraction of the total 
cloud energy at any instant are also shown. 



poses on the hydrodynamic grid helps suppress some of the 
random fluctuations during the interac tion, and also leads t o 
better convergence in resolution tests jPittard et al.|[2009l ). 




Figure 14. As Fig. 1131 but each panel is for models with the 
same shock Mach number: M = 1.5 (top), M = 3 (middle), and 
M = 10 (bottom). The central plot of each panel now directly 
shows the thermal energy of cloud material. Note that the scale 
for the turbulent energy in the top panel differs from that in 

Fig.m 

4.3 Turbulence in shock-cloud interactions 

The interaction of shocks with clouds is likely to be a key 
mechanism for generating turbulence in the ISM, since sub- 
stantial vorticity and velocity dispersion is produced. Tur- 
bulent motions in our simulations can be identified on two 
separate scales: the subgrid turbulence model deals with tur- 
bulent motions and mixing on scales smaller than the size 
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X ^ 10^, with most of the turbulent energy being dissipated 
by f = 10 tec- In contrast, when x ~ 10, the sub-grid turbu- 
lent energy can be roughly constant for durations in excess 
of 15 ice- 



Figure 15. Top: tdrag (for the cloud) and bottom: tmix (for the 
core) as a function of Mach number and x measured from the k-e 
numerical simulations. 



of a grid cell, while larger scale (resolved) turbulent motions 
can be directly measured from the velocity dispersions in 
the axial and radial directions, Sv^ and Svi, respectively. 



4-3.1 The cloud velocity dispersion 

Fig. 1161 shows the time evolution of the velocity dispersions, 
(5uz and Svi. Their behaviours are broadly similar for interac- 
tions with supersonic post-shock flow {M = 3 and M = 10). 
However, the velocity dispersions in the A/ = 1.5 simula- 
tions build more slowly, and are limited to lower maximum 
values, reflecting the weaker interaction. Sv^ peaks at higher 
values with x, due to greater growth of KH and RT insta- 
bilities resulting from the longer drag and mixing timescales 
with X- The ratio Svz/Svr is not strongly dependent on M, in 
contrast to it s variation with x- The latter point was not ed 
previously bv lKlein et al] (|l994l ) and lPittard et all (|2009l ). 



4-3.2 The subgrid turbulent energy 

The time evolution of the energy in sub-grid turbulent mo- 
tions is shown in Fig. 1131 It rises more rapidly and peaks at 
higher values in simulations with high Mach numbers and 
cloud density contrast. When x ~ 10^, it peaks at about 
3 percent in models with M = 3 and M = 10, but peaks 
at only 1 percent when M = 1.5. In models with x = 10 
it peaks at less than 1 per cent, and reaches a maximum 
value of only 0.2 percent when M — 1.5. The peak in the 
sub-grid turbulent energy fraction is relatively narrow when 



4.3.3 Locally averaged velocity dispersions 

The velocity dispersions shown in Fig. [16] are global aver- 
ages over the entire cloud. These can be used to obtain a 
globally averaged energy fraction of the fluctuations. For 
model ml0c3, Svr and Sv^ peak at f « 7.25 tec, with values 
of « 0.09 Wb and ~ 0.23 Ub respectively. H ence, the global 
velocity dispersion, 5v = y/ (2Sv y + 5vi)/3 ~ 0.15 tib. T his 
is about twice the value noted bv lNakamura et al.l (|2006l ) for 
their simulations of radiative clouds with 7 — 1.1. The ratio 
of the energy in (resolved) fluctuations to the kinetic energy 
of the mean flow is (Sv/v)^ ~ 0.2. 

However, such global estimates of Sv and (Sv/v)'^ are 
actually upper limits. Fig. [IT] shows the mean axial veloc- 
ity, Uz, and the velocity dispersions Svz and Sv^, evaluated 
on a more local scale. The tiles in these maps are of size 
0.5 X 0.5 Tc, and contain up to 4096 individual cells from the 
hydrodynamic calculation, over which the displayed quanti- 
ties are averaged. These averages are again mass-weighted 
quantities, but consider material only within each tile - i.e. 
Svz is again calculated using Eq. 23 in Paper I, but {v'^) and 
(tiz) are evaluated only for the material within each tile. 
Note that Svi will be non-zero even when Vr is constant over 
each tile. This is a consequence of the specific averaging for- 
mula. The maps in Fig. 1171 were calculated at t = 5.66 tec 
for clouds with x ~ 10'' hit by a Mach 1.5 (a) and a Mach 
10 (b) shock (for which the corresponding density plots are 
shown in Fig. [2)), and also at f = 7.10 tec for the M — 10 
case (panel c). 

In Fig. I17r a) and (b) , the mean axial velocity is high- 
est at the slip surface between cloud and ambient material, 
and is lowest in the central core of the tail. A velocity gra- 
dient perpendicular to the tail is clearly present. Both the r 
and z— components of the velocity dispersion are also great- 
est at the slip surface. Also shown is logio{5v/v)^ , which 
has a maximum value of —0.1 in Fig. I17f a). and —0.92 in 
Fig. I17r b). Slightly later in the ml0c3 simulation, a large 
fragment breaks away from the cloud core. This fragment is 
responsible for a significant increase in the local and globally 
averaged velocity dispersions, as seen in Figs. I17r c) and 1161 
respectively. The peak local value of {5v/v)^ — 0.6, indi- 
cating that there are regions in the fiow where the swirling 
motions within the gas are almost as fast as its average bulk 
speed. 

Our aim in this section is simply to highlight how 
the velocity dispersion varies on smaller scales in the flow. 
This is of interest given that spatial variations in the ve- 
locity and the velocity dispersion of clouds and tails can 
be p robed using high-spa t ial-resolution observations (see, 
e.g., iMeaburn et"alTll998l : iMeaburn fc BoumisI |2009| ). We 
also wished to draw attention to the fact that while the local 
velocity dispersion may exceed the globally averaged value, 
the latter is in turn an upper limit for the whole cloud. 
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4.4 The mass-loss rate and lifetime 

The time evolution of the core mass and various mass dis- 
tributions have been previously discussed in Section 14.2.31 
Here we note a more detailed examination of this mass-loss. 
In Fig. llSl the rate of mass-loss from the core is compared to 
the analytical formula for hydrodynamic ablation given by 
iHartguist et al.l (| 19861 ). For this comparison we assume that 
the cloud is fully ionized with a radius rc — 2 pc, core density 
Pc = 4 X lO"'^^ g cm~^, and temperature of 8000 K, and is 
in pressure equilibrium with the surrounding gas of density 
4 X 10"^* g cm"^ (i.e. x = 10^) and temperature 8 x 10*^ K. 
Both the cloud and its surroundings are assumed to have so- 
lar abundances (average mass per particle, ^ — 0.61 mn). A 
shock of Mach number 1.5, 3, or 10, then travels through the 
ambient medium at a speed of 650, 1300, or 4300 kms~^ , 

heating the medium to 12, 29, or 260 millio n K. 

The analytical formula proposed by IHartguist et al.l 

l|l986l ) depends on the Mach number of the flow past the 
cloud. If this flow is supersonic (such as when the shock 
Mach number, M > 2.76), then M^b « {McCcf^^ipvY''^, 
while if the Mach number of the flow past the cloud, Mps, 
is subsonic, then there is an additional Mpl^ dependence. 
Mc is the mass of the cloud, Cc is the sound speed within 
it, and p and v are respectively the density and velocity of 
the environment external to the cloud (the post-shock ambi- 
ent /intercloud gas). With the above parameters, Mab ~ 1.1, 
4.9, and 8.2 x 10~^ M0yr~^ for interactions with Mach 
number 1.5, 3, and 10 shocks, and the cloud survives for 
approximately 1.8, 0.41, and 0.24 million yrs. In comparison 
the cloud crushing timescale, tec = 9.5, 4.8, and 1.4 x 10'* yrs, 
so that the cloud survives for about 18 cloud crushing 
timescales before being destroyed (clouds with a lower den- 
sity contrast survive for a smaller multiple of tec)- 

From Fig. [T8]we see that the level of agreement between 
Mab and the numerically determined mass-loss rates from 
the simulations is reasonably good (gaps in the numerical 
curves indicate short periods when the cloud core accretes 
material). Furthermore, the rate of mass-loss from the core 
shows less extreme variations as the Mach number is reduced 
(i.e. the cloud is destroyed in a gentler fashion). However, 
the agreement shown in Fig.[T8]is actually rather fortuitous, 
as we now demonstrate. The lifetime of the cloud, tufc = 
Mc/Mab, can be expressed in units of the cloud crushing 
timescale, tec, as 

1/3 

^lifc / ^cc 



pv 



= 1.6 



pcCc 

pv 



2/3 

1/3 



M 



4/3, 



4/3i 



(13) 



where the term in square brackets should only be used if the 
post-shock flow of density p and velocity v over the cloud is 
subsonic. As M — >■ oo, because of the linear dependence of v 
on M, one finds that tufe/tcc oc M^^^. For constant Pamb and 
Camb (as in the simulations), Eg.lTSlgives tiifn/trr oc (pcCc)^^^, 



-1/2, 



tliic/tcc OC X" 

The ratio of the core lifetime from the numerical mod- 
el^ to the lifetime assuming steady mass-loss at the rate 

^ Defined as the time at which the core disappears - i.e. when ma- 



using the analytical formula of iHartouist" et all (|l986l ) IS 
shown in Fig. 1191 The kink seen in the analytical curves at 
M = 2.76 reflects the switch from Mach number dependent 
mass-loss when the postshock flow is subsonic {Ad < 2.76 - 
i.e. the inclusion of the Mp^^ term in the above equations) 
to Mach number independent mass-loss when the postshock 
flow is supersonic. The displacement of the x = 10^ and 
X = 10'^ curves reflects the x^^^ proportionality noted above. 

The agreement for clouds with density contrasts x ~ 
10'^ hit by shocks with Mach numbers M < 10 is again 
good, as was shown previously in Fig. 1 181 However, there is a 
significant and increasing divergence between the numerical 
and analytical cloud lifetimes as the shock Mach number 
increases past M = 15. For a Mach 40 shock, the numerical 
results indicate that the cloud is "destroyed" by t ~ 8tcc, 
whereas the analytical formula suggests a cloud lifetime of 
20 — 40 tec for X = 10 — 10'', or up to 5 times longer. This is 
because of the M^^^ scaling of Eq.[T3]at high Mach numbers. 
Note that the numerical results are consistent with Mach 
scaling (Section (H)), whereas the analytical formula is not. 
Perhaps even more serious is the disagreement at moderate 
Mach numbers (M < 7) and low density contrasts (x ~ 10). 
In such cases the clou d survives apprec i ably l onger than the 
analytical formula of IHartguist et al.l (| 19861 ) suggests. For 
X = 10 and M = 3, the analytical formula suggests that 
the cloud will survive until t — 4: tec, whereas the numerical 
calculation suggests the cloud actually survives until about 
t = 15 tec (i.e. its mass-loss rate is about 4 times lower than 
the equation in Hartquist et al. would suggest). 

These differences have consequences for previous 
mass-loading calculati ons in which the analyt ical mass-loss 
rate pres cription from Hartguist et al] J 19861) was adopted 
fsee. e.g..[bvson fc Hartguistlll987l:lArthur fc Hennevll 19961: 



Strickland fc Stevens! l2000l : IDv son. Arthur fc HartguistI 
2002h . However, in practice, we believe that the results 



from such calculations will actually change very little, for 
a number of reasons. Firstly, mass-loading causes a flow's 
Mach number to tend towards unity or thereabouts: super- 
sonic flows are slowed as a result of momentum transfer 
and their temperature and sound speed increased through 
frictional heating, while subs onic flows are a c celerat ed by 
mass-loading (see Eq. 5 in IHartguist et al.l (1 19861). and 
the n umerical calculations of Arthur. Dyson fc HartguistI 
(1 19931) ). This means that although use of the formula in 
IHartguist et all 1 19861 ) for a cloud overrun by a high Mach 
number shock or immersed in a high Mach number flow 
results in the initial mass-loading being too slow, the rapid 
decrease in the Mach number of that flow in response to this 
mass-loading means that the analytically determined mass 
injection rate from clouds which are further downstream is 
closer to the correct (numerically determined) value. This 
behaviour is seen, for instance, in the mass-loading of the 
high Mach number pre-termin ation-shock stellar wind ma- 



terial in wind-blown-bubbles (|Pittard. Hartguist fc Dvson 



l2001al ). It is also seen in the simulat io ns of lArthur et al 
1 19931 ) and lArthur. Hennev fc DvsonI (| 19961 ) . where the 
position of the reverse shock of the wind-blown-bubble is 
flxed at the radius of the onset of mass- loading. Secondly, in 



terial originally in the core has at least as much ambient material 
in each computational cell on the hydrodynamic grid. 
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most of these works values for the clouds (such as x, Vc, etc.) 
were not explicitly specified. Instead, the authors simply 
adopt a constant rate of mass-injection into supersonic flow, 
modified by an M"^^^ dependence for subsonic flow. Such 
simulations are therefore "exempt" from the large disparity 
beteen the analytical and numerical results for X ^ 10 ^ind 
M < 7 which we have discovered in this worlQ. 

Having established the existence of a large differ- 
ence between the analytical and numerical lifetimes of 
clouds at high Mach numbers and at low-to-moderate Mach 
numbers for clouds with small density contrasts, we now 
wish to know the underlying cause of this discrepancy. 
For a flow streaming super s onically (i.e. Mpg > 1) past a 
cloud. iHartguist et al.l (| 19861 ) assumed that mass- loss occurs 
largely as a result of low pressure regions on the surface of 
the cloud. The mass-loss rate and lifetime of the cloud (in 
cgs units) is then found to be Mach number independent. 
However, this leads to a Mach number dependent lifetime for 
the cloud when expressed in units of tec (Eg. I13p . since the 
dependence is introduced by tec. In contrast, the numerical 
simulations presented in this work (and earlier simulations 
in the literature) show that, for clouds with x ^ 10^, mate- 
rial is initially stripped from the sides of the cloud through 
the development of KH instabilities, while the core remains 
relatively intact. At later times, and at early times for clouds 
with X ^ 10, the instabilities break up the cloud into sev- 
eral large fragments. It seems reasonable to conclude, there- 
fore, t hat the mass-loss process identified bv lHartguist et al] 
(|l986l ) is in fact not the dominant mode of mass-loss from 
deformable clouds subject to large velocity shear. 

In contrast, Eq. [T^] shows that the destruction time 
due to long wavelength KH instabilities has only a weak 
dependence on the shock Mach number. One finds that 
fKH = 2.4 tec when M = 1.5, and fxH = 1.33 tec when 
M — 40. These estimates are roughly 6 times shorter than 
the numerically determined lifetimes at high shock Mach 
number (Af > 7). The scaling of the cloud lifetime (in units 
of tec) with M is also weaker than is observed from the 
simulations (see Fig. I18|) . For instance, it predicts that the 
lifetime of a cloud hit by a Mach 1.5 shock is about 1.8 times 
longer than in a Mach 40 interaction. Nonetheless, the in- 
sensitivity of Eq. [T2]to the cloud density contrast, x, is con- 
sistent with the numerical results, where only a weak (and 
non-montonic) dependence is seen (Fig. 118(1 . Given these 
observations, it seems reasonable to associate the mode of 
cloud destruction with KH instabilities, but it is clear that 
the lifetime of the cloud, tutc ~ 6tKH. This expression is 
within a factor of 2 or so of the numerically determined 
lifetime of clouds with 10 < x < 10^ hit by shocks with 



* Obviously other works examining different types of mass- 
loading in diffus e sources, for example through thermal 
cond u ction (e.g. [M cKcc fc Ostrikcr ' l977l: IC hieze &: Lazarefj 
198ll; IWhite fc Lona 1991: Pittard Hartguist fc DvsonI l200ld : 
Pittard. Hartquist fc Ashmorell2003ai:|Pittard et al.ll2003hj. 120041) 
or photoevaporation (e.g . iGarci^Airedondo'^rthur'irifa 
I2OO2I: iPittard et al] 120051) . or where the mass-lo ading rate is 
unrelated to the properties of the local flow (e.g. [ S mithl [l 99^: 
Dyson. Williams fc Redmanlll995l: IW illiams. Hartguist fc DvsonI 
I995I: IWilliams. Dyson fc HartquistI Il999l). or is not 
from embedded clouds (e.g . IWeayer et al.l Il977l : 
iToniazzo. Hartquist fc Durisenll200ll) , remain unaffected. 
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Figure 17. Snapshots of (t^z), ^v^, Sv^, and logio{Sv/v)'^ . In each 
pixel in these maps the mean axial velocity, the velocity disper- 
sions, and the ratio of the energy in (resolved) turbulent motions 
to the kinetic energy of the mean flow are calculated from sum- 
ming over many smaller hydrodynamic cells. The maps were con- 
structed at t = 5.66 tec from simulations with x = 10'^ s-i^d ^) 
M = 1.5 and b) M = 10. Panel c) shows the M = 10 simula- 
tion at t = 7. 10 tec I which is near the time of the peak value of 
the globally averaged (Siiz (see Fig. I16f ). The calculations were 
only performed in regions containing cloud material (i.e. where 
the colour in the above plots is other than white). 
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Figure 16. Top panels: Time evolution of the eloud veloeity dispersion in the radial direction, cloud i for various Mach numbers and 
density contrasts: (a) x = 10> (b) X = W^, (c) x = W^. Bottom panels: as the top panels, but for the velocity dispersion in the axial 
direction, cloud- 



Table 4. The parameters from least squares fits to tdrag/*cc for 
the cloud (see Fig. 115b ). The final column lists the sum of the 
residuals. In each case the fitted function has the form tdrag/^cc = 
aM' + c. Each fit has 4 degrees of freedom. 



Table 6. The parameters from least squares fits to tnfe/tcc for 
the cloud. The final column lists the sum of the residuals. In each 
case the fitted function has the form iufc/tcc = aM^ + c. Each fit 
has 4 degrees of freedom. 



S res2 



10 5.11 ± 0.35 -2.33 ±0.16 0.97 ± 0.03 0.0487 
10^ 7.94 ±0.84 -2.46 ±0.25 3.08 ±0.07 0.043 
10^ 5.67 ±0.86 -0.66 ±0.34 4.84 ± 1.04 0.956 



10 25.2 ±10.1 
10^ 20.3 ±4.94 
10=* 66.6 ±21.6 



S res^ 



-1.89 ±1.06 10.1 ±2.00 13.9 
-1.62 ±0.56 8.26 ± 1.06 6.52 
-2.82 ±0.76 10.2 ±1.20 16.3 



Table 5. The parameters from least squares fits to tmix/icc for 
the cloud (see Fig. 115b ). The final column lists the sum of the 
residuals. In each case the fitted function has the form tmix/*cc = 
aM^ + c. Each fit has 4 degrees of freedom. 



X 


a 


b 


c 


E res^ 


10 


13.9 ± 1.67 


-2.90 ±0.28 


6.00 ±0.09 


0.086 


102 


6.93 ± 0.69 


-1.59 ±0.23 


4.67 ±0.15 


0.135 


10^ 


10.7 ± 1.95 


-1.10 ±0.43 


6.11 ±0.96 


2.81 



M > 1.5, with the biggest discrepancy at low Mach num- 
bers. The parameters of the fits to the numerical results are 
noted in Table [6l 



5 SUMMARY AND CONCLUSIONS 

This is the second of a series of papers investi gating the tur- 
bulen t destruction of clouds. In our first paper jPittard et al.l 
l2009l ) the focus was primarily on the much faster evolution 
and dispersal of a cloud over-run by a shock with a highly 
turbulent post-shock flow. Here we have performed a de- 



tailed examination of the Mach number dependence of the 
destruction of a cloud by an adiabatic shock. We have used 
a hydrodynamical code which incorporates a sub-grid com- 
pressible k-e turbulence model in an attempt to calculate 
the properties of the turbulence and the resulting increase 
in the transport coefficients. 

We find that the most significant differences in the na- 
ture of the interaction occur between those where the am- 
bient post-shock flow is subsonic with respect to the cloud, 
and those where it is supersonic (the latter requires that 
the shock Mach number Af > 2.76). For this reason, the 
interaction of a Mach 3 shock is more akin to a Mach 10 
shock than a Mach 1.5 shock. At high Mach numbers the 
post-shock conditions are virtually independent of the Mach 
number and the so-called "Mach scaling" is obtained. Mach 
scaling appears to hold also when using the k-e turbulence 
model. 

For weak shocks (M < 2.76) the interaction is much 
milder and the following differences are observed: i) the 
postshock flow is subsonic with respect to the cloud, so a 
bowwave rather than a bowshock forms ahead of the cloud, 
ii) the compression of the cloud is more isotropic, iii) a 
weaker vortex ring is produced, iv) the smaller velocity dif- 
ference at the slip surface around the cloud limits the KH 
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Figure 18. Mass loss rate for models with x = 10'' ^nd a) M = 
1. 5, b) M = 3, and c) M = 10, compared to the mass-loss formula 
of iHartQuist et al. I j 19861) . Gaps in the curves indicate periods 
where the cloud core temporarily accretes material. 



and RT instabilities and reduces the peak turbulent energy 
fraction of the flow, v) it takes much longer for the cloud to 
be mixed into the surrounding flow and for it to accelerate to 
the intercloud postshock speed, and vi) mass stripped from 
the cloud does not as readily form a long tail (the set-up 
time is longer). We further flnd that a prominent tail only 
forms if X > 10^ 

Our m ost important finding is that the analytical pre- 
scription in lHartguist et ahl (|l986h for the ablative mass-loss 
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Figure 19. Comparison of the cloud lifetime from the numerical 
simulations (defined as the time when the "core" disappears) and 
the lif etime obtained from the analytical result of lHartouist et aU 
lll986l) . as a function of the shock Mach number, M, and the cloud 
density contrast, x- 

rate of a cloud in an external flow predicts cloud lifetimes 
which are in disagreement with numerically determined val- 
ues. For instance, the predicted lifetimes are a factor of 
2 — 5 times too long for clouds with 10 < x < 10'^ hit by a 
Mach 40 shock, while they are about 4 times too short for 
clouds with X = 10 hit by an M = 3 shock. The reason for 
these discrepancies ap pears to be due to the assumption in 
iHartguist et al.l(| 19861 ) that the mass-loss is mostly driven by 
pressure gradients around the cloud. Instead, we show that 
the cloud lifetime is more closely related to the timescale 
for large scale KH instabilities, though it is about 6 times 
longer than the latter. We argue, however, that the rapid re- 
duction in the Mach number of hypersonic flows subject to 
mass-loading means that previous work in the literature is 
unlikely to greatly change if repeated using a more accurate 
mass-loss rate prescription. 

In future work we will extend our investigation to three 
dimensions, examine the interaction of a dense shell with 
a cloud, and will compare synthetic signatures of the in- 
teraction to the types of diffuse sources mentioned in the 
introduction. 
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